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Abstract 

A simplified single rate expression for hydrogen combustion and nitrogen oxide production 
was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete 
chemical mechanism over the entire operating space. These times are then correlated to the 
reactor conditions using an exponential fit. Simple first order reaction expressions are then used 
to find the conversion in the reactor. The method uses a two-time step kinetic scheme. The first 
time averaged step is used at the initial times with smaller water concentrations. This gives 
the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, 
and pressure. The second instantaneous step is used at higher water concentrations 

o rv 

(>1x10 moles/cc) in the mixture which gives the chemical kinetic time as a function of the 
instantaneous fuel and water mole concentrations, pressure and temperature (T4). The simple 
correlations are then compared to the turbulent mixing times to determine the limiting properties 
of the reaction. 

The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for 
each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate 
the necessary chemical kinetic times. This time is regressed over the complete initial conditions 
using the Excel regression routine. Chemical kinetic time equations for H 2 and NO x are obtained 
for EE/air fuel and for the H 2 /O 2 . 

A similar correlation is also developed using data from NASA’s Chemical Equilibrium 
Applications (CEA) code to determine the equilibrium temperature (T4) as a function of overall 
fuel/air ratio, pressure and initial temperature (T3). High values of the regression coefficient R 
are obtained. 


Introduction 

Hydrogen is the fuel of the future. Iceland is converting their economy to hydrogen and 
several countries (United States, Germany, Australia, Mexico, China, Canada, Japan, and Russia) 
are strongly interested in the technology because of global warming concerns from CO 2 and the 
unlimited supply of hydrogen from water conversion. It is important to have a simple, quick 
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method of computing the chemical rates for prediction of times of conversion and for the design 
of combustors. Accurate chemical kinetic mechanisms exist for hydrogen reaction, but they 
involve as many as 2 1 chemical steps with the production of intermediate species like OH, HO 2 , 
H and O atoms with H 2 O 2 as an intermediary. We have used the H 2 mechanism of S.L. Hwang 
et al. (ref. 1), as the basis of this report. No changes in the hydrogen mechanism were made in 
switching from H 2 /Air to H 2 /O 2 . We have taken a forty-two step mechanism for the production 
of nitrogen oxides from Gri-Mech 2.1 without the carbon reactions, used in reference 2, and 
correlated the formation rate of nitrogen oxides in the presence of hydrogen by a similar two- 
time step method. Any detailed mechanism could have been used to compute the chemical 
kinetic time constants. 

Large mechanisms with many intermediate species and very fast radical reactions cause the 
equations to be stiff (extremely fast compared to the overall rate, requiring a large number of 
small time steps), making them very difficult to integrate. Calculations for these extensive 
mechanisms are repetitive and complex. Using the simplified kinetic scheme developed here to 
calculate the two chemical kinetic times greatly reduces the amount of time required to compare 
kinetic reaction times with turbulent mixing times and will reduce the time required to obtain a 
converged solution. The advantage of extracting the chemical kinetic time for only the species of 
interest from a detailed computation is that we have only the differential equations of interest to 
solve, resulting in a much smaller set of equations. 

The hydrogen reactions are about seven times as fast as hydrocarbon fuels and the reactions 
require small time steps to integrate accurately. We have used the NASA Glenn GLSENS 
kinetics code (ref. 3) to integrate the progress of the reaction accurately. The classical methods 
were used to compute over two thousand cases at constant T and P. Then the results were 
correlated with a simple two-time step approach to predict the chemical kinetic time constants. 
The derived method involves only the principal species H 2 , and H 2 O and NO x for nitrogen 
oxides, species which are predicted routinely by computer codes. We do not require tracking the 
intermediate free radicals, as the kinetic times in step one before the radicals become large are 
practically constant and in step two are related to the instantaneous H 2 , O 2 , and H 2 0 
concentrations . 

We have also curve fitted the equilibrium temperature T4 using data generated by the NASA 
Chemical Equilibrium Application code (CEA) to an exponential equation. 

The NASA detailed kinetics code (GLSENS) (ref. 3), was used to integrate the system of 
equations at constant temperature and pressure, at over 2000 initial conditions to derive the rate 
expressions. We have massively correlated the output from GLSENS, into simple exponential 
expressions for the chemical kinetic times using the Excel regression routines. It may be 
reasoned that the presented equations are only as good as the overall mechanism that calculates 
the data. However, performing the calculations in the conventional manner is also only as good 
as the mechanism equations and constants that go into them. 

The two-time step approach was used successfully by Molnar and Marek in reference 2 for 
Jet-A and methane. The two-time step approach is used because the free radicals OH, etc., 
catalyze the reactions and act as important intermediates. The principle reaction for water 
production is OH through the reaction 

H2 + OH = H20 + H (1) 
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We think that the more water that is present in the reactor then the more OH would be 
present accelerating all of the reactions. So H 2 0 can be used as an indicator of the OH present. 

The report will present the section Mixing/Kinetics, Model Equations, Equilibrium Fits, 
Determination of the Chemical Kinetic Time, and the Chemical kinetic time equations for both 
the H 2 /air and the H 2 /0 2 mixtures. The principal difference in the conditions is in the maximum 
range of pressures and temperatures considered, see Table Input, and the lack of the diluent 
nitrogen for the H 2 /0 2 condition. Comparisons of the correlations for different fuels and a 
comparison of NO x production data from reference 2 is given. 


Mixing Times Versus Chemical Kinetics Times 

Reaction rates are kinetically limited at low temperatures and mixing limited at very high 
temperatures. According to the Magnussen model (ref. 4), the fuel oxidation rate will be 
determined by the maximum of either the chemical kinetic time or the turbulent mixing times of 
the fuel and air. However, for large numerical solutions it is very tedious to use detailed classical 
calculations to compare both the kinetic and turbulent mixing times to determine the limits of the 
reaction. Detailed chemical kinetic schemes are extremely time consuming for two- and three- 
dimensional computer calculations for combustors. 

The Magnussen model (ref. 4) proposes that the maximum of either the turbulent mixing or 
the chemical kinetic times will be the limiting factor of a chemical reaction. This model could be 
explored by numerically calculating both times to compare them. However, using detailed 
mechanisms, this is a long and tedious process. The calculations would be extremely 
complicated for the detailed chemical kinetic time. By using the equations presented here to 
determine the chemical kinetic times and using conventional numerical methods to determine 
mixing times, the Magnussen model can be applied in a much more convenient way. 

Net rate m r = min (^-y fud , A£ ^ ox>,gew ,m kinetic ) (2) 

k k rj- 


Where — equals the turbulent mixing time, r m , with k being the turbulent kinetic energy, s 

As 

is the dissipation rate, y is the mass fraction, and ry, is the stoichiometric coefficient written on a 

y fuel 

mass fraction basis. The mixing constant, A, is usually given as 4.0. The factor — — - is the 

^kinetic 

chemical kinetic time x c computed in this report from the correlations presented. 

In order to obtain the chemical source termer,. , a comparison is made of the mixing rate, 

— and the chemical kinetic rate — , and the lowest rate or the longest time is used in the 
r r 

m y c 

expression (see fig. 1). This may also be represented by the following relationship: 

t = max ( T m , t c ) (3) 
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With the approach derived here, a simple direct comparison can be made between the mixing 
and chemical kinetic times and the minimum rate used for the computation. 


Model Equations 

The following equations can be used to model the chemical system. 


H 2 + 1/2 0 2 


^Fuel 

► h 2 o 


(SI) 


zno 

14 N 2 + '/ 2 0 2 ► NO 


(S2) 


The following first order reaction was used to represent the rate of fuel burning. We compute 
the net reaction time so we do not worry about the forward and reverse processes, only the net 
process. We assume that the fuel goes to zero at the end of the reaction and that the NOx 
increases linearly and it is far from equilibrium. (In this report, t and z are given in milliseconds, 

TYIS * CC 

except Tno x which is in , while concentrations are given in gmoles/cc): 

gmol 


dm 

dt 


cm 

t H 2 


(4); SO 


cHl 


t H2 ~ 

I 


dm ' 

dt j 


OR if Th 2 is constant: 


Hi = Hie 


-t 
Th 2 


For a constant zm, the fuel concentration is then represented by a simple exponential 
decay expression, where H 2 o is the initial fuel concentration. 

The nitrogen oxide formation rate, a species important for combustor emissions, was 
modeled as a simple zero order expression. 


dNOx 1 t 

— ~T — = j (5) OR if Tnox is constant NOx = (6) 

Clt X /r 

NO X L NOx 


* , . r MS -CC 

t nox has units ot 

gmol 
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Determination of Chemical Kinetic Time 


Using equations (4) and (5) above, the dependent chemical kinetic times (tfuci , tnox) can be 
easily computed and correlated against several independent variables using the Excel regression 
routines. The integration was performed for 1256 cases for the H 2 /air and 2808 cases for the 
H 2 O 2 conditions shown below: 

TABLE-INPUT: INITIAL CONDITIONS FOR COMPUTATIONS 

H2/Air Inlet Conditions: 

Equilibrium H2/Air 

H2(g) = 294.44 K 

Air =294.44 K, 589.89 K, 922.22 K, 1255.55 K 
Pressure = 1,20, 40, 60 atm 

Lean equivalence ratios = 0.2 to 1.0 increments of 0.1 
Rich equivalence ratios = 1.05, 1.1 to 2.0 increments of 0.1. 

Initial conditions for chemical kinetic times for H2/Air 
Temperature = 1000., 1500., 2000., 2500. K 
Pressure = 1, 10, 20, 30, 40 atm 
Lean equivalence ratios = 0.3 to 1.0 increments of 0.1 
Rich equivalence ratios = 1.05, 1.1 to 2.0 increments of 0.1 

H2/02 Inlet Conditions: 

Equilibrium H2/02 

H2(g) = 30., 100., 200., to 500. K increments of 100. K 

02(g) = 50., 100., 150. K 

Pressure = 1, 50, 100, 200, 300, 400 atm* 

Lean equivalence ratios = 0.1 to 1.0 increments of 0.1 
Rich equivalence ratios = 1.05, 1.1, 1.2 to 3.0 increments of 0.2 
Initial conditions for chemical kinetic times for H2/02 

Temperature = 1000., to 4000., increments of 500, then 5000. K 
Pressure = 0.5, 1., 10., 50., 100., 200., 300., 400. atm 
Lean equivalence ratio = 0.1 to 1.0 increments of 0.1 
Rich equivalence ratio = 1.05, 1.1, 1.2 to 3.0 increments of 0.2 
* These conditions are for advanced rocket combustors (ref. 5). 

on 

Time steps for step 1 was integrated to a water concentration of 1x10 and time 
step 2, the instantaneous chemical kinetic times were computed to 0. 1 ms. The output for 
time step 2 is every integration step after the water concentration was above 1x10 20 , with 
up to 20 points tabulated per case. The total number of points could exceed 20 000. 
Microsoft Excel (Microsoft Corporation, Redmond, WA) can handle up to 65 000 points. 

The supplied kinetic mechanism was used for the full range of conditions without 
modification. Even though the mechanism was tested to a temperature of 3 000 K and 4 atm, we 
felt that the mechanism would extrapolate and give reasonable results for design since nothing 
else was available at the high temperature and high pressure conditions. The differences between 
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the th/air and the H 2 /O 2 is the absence of nitrogen in the H 2 /O 2 set, and in H 2 /O 2 the maximum 
temperatures were twice that of H 2 /air and the maximum pressures were an order of magnitude 
higher. We computed and correlated the chemical kinetic times over the full range of conditions. 

For all time steps, calculations were performed isothermally using GLSENS for each 
condition over a time of 0 to 3 ms. By computing the progress isothermally, the chemical rate 
constants were nearly fixed and averaged over the total time for step one and instantaneous 
chemical kinetic times were computed for step two: 

—!— = — [ (7 

Taverage ttotal T instantaneous 


The chemical kinetic time was determined as a unique value of temperature, pressure, and 
instantaneous mole fractions of fuel and water. GLSENS computes the cumulative rate of 
reaction for each species from all equations in the mechanism, so it is a simple matter to then 
compute the chemical kinetic time for each species. 

dcH2 _ v dcH2 1 
dt ' total ~ dt I i 


For the fuel equation (4) the chemical kinetic time is given as 


r 


H 2 


cHi 


Z 


dHi' 

dt j 


( 9 ) 


This simple calculation was done using additional steps in Subroutine Out2 in the GLSENS 
code (see appendix B). Values for the chemical kinetic time were calculated for each 
concentration at each output time and each set of conditions. For time step 1, the trapezoidal rule 
(using 1/r) was used to average the chemical kinetic time to calculate the best value for each set 
of conditions and the final numbers regressed over the complete set of cases to obtain the final 
correlation. 

A correlation could then be developed that determines the chemical kinetic time as a 
function of the initial overall cell fuel/air ratio, pressure and temperature. The data was correlated 
using the same method as previously mentioned for the equilibrium equations. Two correlations 
for each step for each of the two species, one for the lean side and one for the rich side, were 
obtained. This results in a total of 4 H 2 correlations for F^/air and 2 H 2 /O 2 correlations. As OH is 
formed, the reaction rates get faster and the chemical kinetic time gets smaller. The chemical 
kinetic time was correlated using two different steps to increase the accuracy of the calculation. 
Step one is an average chemical kinetic time taken over 3.0 ms. Step two is an instantaneous 
value that depends on the instantaneous amounts of fuel, water, and oxygen. Since we were not 
tracking the radical concentrations, we were using H 2 O to indicate the state of the radical 
species. In other words for the reactions 

H2 + OH = H20 + H (10) 
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N + OH =NO + H etc. 


( 11 ) 


The radicals could be correlated with H 2 O because of the first reaction. 

The user would switch from step one to step two when the molar concentration of water is 

—20 

greater than 1^10 moles/cc. This value was chosen arbitrarily as a simple switch to transition 
when the rates accelerated. Step two can not be used with small concentrations of water because 
if the value of the water concentration was zero, the entire correlation time would go to zero and 
the rate would be computed as infinity. 


H 2 Kinetic Scheme For Step 1 (Average) and 
Step 2 (Instantaneous) Methods 

The following is GLSENS input for the 21 step, 8 species mechanism from reference 1 and 
the 42 step NO equations from GRI-Mech 2.1, without the carbon reactions from reference 2 that 
were used for the H 2 calculations for a total of 63 mechanism steps and 15 species. 

TAPE 

NEW H2 MECHANISM OF S.M.Hwang, et. al. (reference 1) 

&RTYPE GLOBAL=.TRUE.,GRONLY=.FALSE. &END 


H 

02 

— 

OH 

0 

6.73E+15 

-0.499 

26671.6 

0 

H2 

= 

OH 

H 

5.06E+04 

2.67 

6289.2 

OH 

H2 

= 

H20 

H 

2.16E+08 

1.510 

3430.0 

0 

H20 

= 

2.00 

H 

4.51E+04 

2.7 

14551.4 

0 

O 

= 

02 

M 

1.00E+17 

- 1 . 

0.0 

THIRDBODY 







H20 

18.5 

H2 

2.9 

02 

1.2 

END 


H 

H 

= 

H2 

M 

6.4E+17 

- 1 . 

0.0 

THIRDBODY 







H20 

12.0 

H2 

4.0 

H 

26. 

END 


H 

0 

- 

OH 

M 

6.2E+16 

-0.6 

0.0 

THIRDBODY 







H20 

5. 

END 






H 

OH 

- 

H20 

M 

8.4E+21 

-2. 

0.0 

THIRDBODY 







H20 

16.25 

H2 

2.5 

END 




H 

02 

= 

H02 

M 

5.55E+18 

-1.15 

0.0 

THIRDBODY 







H20 

21.3 

H2 

3.33 

02 

1.33 

END 


H02 

H 

- 

OH 

OH 

8.4E+13 

0.0 

635.0 

H02 

H 

- 

H2 

02 

2.5E+13 

0.0 

693.1 

H02 

H 

- 

H20 

0 

5.0E+12 

0.0 

1410.1 

H02 

0 

- 

02 

OH 

2.0E+13 

0 . 

0 . 

H02 

OH 

= 

H20 

02 

2.0E+13 

0 . 

0 . 

H02 

H02 

- 

H202 

02 

1.3E+11 

0 . 

-1630.0 

H02 

H02 

= 

H202 

02 

4.2E+14 

0 . 

11999.5 

M 

H202 

= 

2.0OH 

M 

1.2E+17 

0 . 

45506.7 

THIRDBODY 







H20 

6.0 

02 

0.78 

END 




H202 

H 

= 

H02 

H2 

1.7E+12 

0 . 

3776.3 
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H202 

H 

- 

H20 

OH 

1.0E+13 

0. 

3585.1 

H202 

0 

- 

H02 

OH 

2.8E+13 

0. 

6405.4 

H202 

OH 

= 

H20 

H02 

1.0E+12 

0. 

0. 

H202 

OH 

= 

H20 

H02 

5.75E+12 

0. 

9559.8 

N 

NO 

= 

N2 

O 

3.500E+13 

.000 

330.00 

N 

02 

= 

NO 

O 

2.650E+12 

.000 

6400.00 

N 

OH 

= 

NO 

H 

7.333E+13 

.000 

1120.00 

N20 

0 

= 

N2 

02 

1.400E+12 

.000 

10810.00 

N20 

0 

= 

2.0 

NO 

2.900E+13 

.000 

23150.00 

N20 

H 

= 

N2 

OH 

4.400E+14 

.000 

18880.00 

N20 

OH 

= 

N2 

H02 

2.000E+12 

.000 

21060.00 

M 

N20 

= 

N2 

O 

1.300E+1 1 

.000 

59620.00 

THIRDBODY 







H2 

2.00 

H20 

6.00 

N2 

0.625 

END 


H02 

NO 

= 

N02 

OH 

2.1 10E+12 

.000 

-480.00 

NO 

O 

= 

N02 

M 

1.060E+20 

-1.410 

.00 

THIRDBODY 







H2 

2.00 

H20 

6.00 

N2 

0.625 

END 


N02 

O 

= 

NO 

02 

3.900E+12 

.000 

-240.00 

N02 

H 

- 

NO 

OH 

1.320E+14 

.000 

360.00 

NH 

O 

= 

NO 

H 

5.000E+13 

.000 

.00 

NH 

H 

- 

N 

H2 

3.200E+13 

.000 

330.00 

NH 

OH 

- 

HNO 

H 

2.000E+13 

.000 

.00 

NH 

OH 

- 

N 

H20 

2.000E+09 

1.200 

.00 

NH 

02 

- 

HNO 

O 

4.610E+05 

2.000 

6500.00 

NH 

02 

= 

NO 

OH 

1.280E+06 

1.500 

100.00 

NH 

N 

= 

N2 

H 

1.500E+13 

.000 

.00 

NH 

H20 

= 

HNO 

H2 

2.000E+13 

.000 

13850.00 

NH 

NO 

= 

N2 

OH 

2.160E+13 

-.230 

.00 

NH 

NO 

= 

N20 

H 

4.160E+14 

-.450 

.00 

NH2 

O 

= 

OH 

NH 

7.000E+12 

.000 

.00 

NH2 

O 

= 

H 

HNO 

4.600E+13 

.000 

.00 

NH2 

H 

= 

NH 

H2 

4.000E+13 

.000 

3650.00 

NH2 

OH 

= 

NH 

H20 

9.000E+07 

1.500 

-460.00 

H 

NO 

= 

HNO 

M 

8.950E+19 

-1.320 

740.00 

THIRDBODY 







H2 

2.00 

H20 

6.00 

N2 

0.7 

END 


HNO 

O 

= 

NO 

OH 

2.500E+13 

.000 

.00 

HNO 

H 

= 

H2 

NO 

4.500E+1 1 

.720 

660.00 

HNO 

OH 

= 

NO 

H20 

1.300E+07 

1.900 

-950.00 

HNO 

02 

= 

H02 

NO 

1.000E+13 

.000 

13000.00 

NH3 

H 

= 

NH2 

H2 

5.400E+05 

2.400 

9915.00 

NH3 

OH 

= 

NH2 

H20 

5.000E+07 

1.600 

955.00 

NH3 

O 

- 

NH2 

OH 

9.400E+06 

1.940 

6460.00 

NH + 

NO 

= 

N2 

+ OH 

2.00E+15 

-.8 

0. 

H 

02 


H20 

= H02 

H20 



1. 

1. 

1. 



9.380E+18 

-.760 

.00 

H 

02 


N2 

= H02 

N2 



1 . 

1 . 

1 . 



3.750E+20 

-1.720 

.00 


The hydrogen mechanism was not changed for the H2/O2 calculation. Note, some reactions 
are bimolecular and some are trimolecular expressions so the kinetics code GLSENS was used 
instead of LSENS. However, the GLSENS code follows the method of LSENS developed by 
Radhakrishnan (ref. 6 ). 
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H 2 /Air Equilibrium 


Equilibrium correlations were generated by using Microsoft Excel to perform a multivariate 
linear regression on the large data set generated by the CEA program of reference 7. (The 
detailed procedure describing the regression used for both equilibrium and finite rate chemical 
times can be found in appendix A). CEA has a plot option for direct tabulation of the output data, 
for ‘f/a, P, T4 H20, NO’. Although the equilibrium correlations were not usually used in the 
calculation of the chemical kinetic times, we feel that these equations could still be very useful 
for other calculations of T4 in well stirred reactors or plug flow calculations. Table 1 shows the 
H 2 equilibrium correlations for T4, and NO x for both the lean and rich cases. The exponential 
function did not correlate the complete data over the total equivalence ratio from 0.2 to 2.0, 
because of the peak in the values at the stoichiometric conditions. We chose to correlate the 
values from lean = 0.2 to 1 .) and then rich from ((|) 0 >1 .0 to 2.0). A similar approach was used 
for the chemical kinetic times. Figure 2 is a parity plot showing the strength of the lean T4-H 2 /air 
equilibrium correlation. This plot shows a minimal amount of scatter, mostly at an equivalence 
ratio of 1.0, indicating a good correlation (R values greater than 0.9). R is based on the 
logarithm of the variables (see appendix A). Using the logarithms, the errors are relative errors 
so errors at the small times are weighted equally to errors at large times. An R 2 of one is a perfect 
correlation. Although the scatter looks large for some parity plots, the mean error looks good 
over the 1000 points. The correlation is very good over a variable range of seven orders of 
magnitude. Note that the units of NO x are mole fraction. This parity plot is typical of all the 

2 ^ 

variables because of the good R values obtained. The lean T4 and NO x equations were 
correlated over the complete temperature, pressure and phi range of independent variables in the 
Table Input. 


T 4 = A(P) a (T 3)\(/)o) c 


( 12 ) 


NOx eq (mole fraction) = A(P) a (<fb) b exp 


E_ 

74 


(13) 


Although these equations seem complex, they are simple to program and use. 
Calculations can even be done on a hand calculator. 


Table 1.— EQUILIBRIUM Lean and Rich H 2 /Air Correlations 
T4 eq Equations 



A 

Pressure 

T3 

<k> 

Rsquared 

Lean T4 

476.9 

0.00419 

0.261 

0.432 

0.936 

Rich T4 

1217.0 

0.00727 

0.120 

-0.243 

0.962 


No x eq Equations 



A 

Pressure 

4*0 

1/T4 

Rsquared 

Lean NO x eq 

0.616 

-0.0153 

-0.810 

-11382.3 

0.977 

Rich NO x eq 

21168. 

-0.475 

-4.080 

-39583.4 

0.991 


(14) 

(15) 


(16) 

(17) 
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Chemical Kinetic Times for H 2 /Air 


Step One Equations for H 2 /Air 

The following form of equation was used for the H 2 and NO x step one correlations: 

T = A(PY(,^) b tx V j (18) 

where r is the chemical kinetic time in milliseconds, P is pressure in atm, f/a is the initial or 
overall mass fuel air ratio as in reference 2, and T is the temperature in Kelvin. The correlation is 

on 

switched to Step Two when the water molar concentration is greater than 1x10 moles/cc, see 
the next section. 

The coefficients for each of the parameters in the correlations may be found in Table 3. 
R-square values have been included to demonstrate the strength of the correlation. R-squared is a 
measure of the error in the model; an R-squared value of one is ideal. 

Since the radical and intermediate species concentrations are very small, the initial or 
overall fuel air ratio can be determined from the molar concentration of cfT, cH 2 0 and 
cN 2 moles/cc as follows: 

f la = (cHi + cHiO) *2* — ; and (|)o=f/a / 0.0292 (19) 

29 * cN 2 

Table 2. — Step One H 2 /Air Chemical Kinetic Time Correlations 

A a E b Rsquared 

H 2 /Air Lean" 4.671 e!T~ -1.036 25068~ 0.539 0.994~ 

H 2 /Air Rich 3.223E-5 -1.032 24223. 1.832 0.951 

No x Lean 2.36E5 -1.670 45085. 0.472 0.647 

No x | Rich | 1.831 | 1.484 | 49802. | -5.91 | 0.397 

All of the results presented here are correlated over the complete conditions listed in the 
Table Input, initial f/a, temperature, and pressure of the reactor. Parity plots for the lean step one 
correlations have been created and may be found in figures 4 to 5. The x-axis contains values for 
the chemical kinetic time generated by the full mechanism GLSENS at each condition. The 
y-axis contains values calculated using the chemical kinetic time correlations above at the same 
set of conditions. This parity plot demonstrates how close the calculated value is to the full 
mechanism value and is a good measurement of the strength of the correlations. The rich parity 
plots are similar to the lean plots shown here and are not repeated. 

Step Two Equations for H 2 /Air 

The following form of equation was used for both the Tau H 2 step two fuel and NO x 
correlations except for different constants. This form of the equation produced the best fit: 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 
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(24) 


t = A(P) a (cH2) b (cH20) c (<f>o) d expy 


where P is pressure in atm, cH 2 is the instantaneous molar concentration of cH 2 and cH 2 0 is 
the instantaneous water concentration. We have preceded the symbol with a “c” to indicate molar 
concentration moles/cc was used and a “y” for mole fraction. Molar concentration is usually a 
primary variable of combustor computer codes. T is the mixture temperature in Kelvin. We have 
correlated to only the major species hoping that cH 2 0 will track the minor species (OH, H, O, 
etc.) to allow good overall correlation and easy to use equations. 

The coefficients for each parameter are given in table 3. Parity plots for the step two 
lean H 2 correlations can be found in figures 6 and 7. These figures show considerable scatter for 
the fuel and small amount for NO x . The hydrogen mechanism used is difficult to model with the 
current equations as indicated by the scatter. However the correlations are good for design and 
large numerical calculations where kinetics are important. 


Table 3 

I. — Constants in Equation (24) Step 1 

rwo H 2 /Air Kinetic Time Correlations 



A 

a 

E 

b 

c 

d 

Rsquared 

H 2 /Air 

Lean 

1.89E-14 

0.605 

18497. 

-0.803 

-0.630 

2.41 

0.604 

H 2 /Air 

Rich 

5.91 E-34 

3.80 

26071 . 

-3.80 

-0.511 

8.24 

0.431 

NO x 

Lean 

1.45E-10 

-0.0092 

48747. 

-0.479 

01 .042 

2.052 

0.967 

NO x 

Rich 

1.26E-19 

1.566 

49662. 

-2.01 

-0.932 

5.99 

0.869 


(25) 

(26) 

(27) 

(28) 


H2/O2 Equilibrium Correlations 

The H 2 /0 2 presentation closely follows the H 2 /air presentation. The lean and rich H 2 /0 2 
equilibrium correlations can be found in table 4, respectively. A parity plot for the lean H 2 /0 2 T4 
equilibrium correlation can be found in figure 8. T3 is the initial temperature for H 2 and 0 2 . 
whose state were taken to be gaseous. These temperatures for rockets is considerably lower than 
for H 2 /air combustors. 


T 4 = A(P) a (0 o ) b (T3H2) c (T 3 o2) d 


(29) 


Table 4.— EQUILIBRIUM H 2 /Q 2 Correlations 



A 

a 

b 

c 

d 

Rsquared 

Lean T4 

1.46E-109 

0.0248 

0.393 

8.96 

58.4 

0.91 

Rich T4 

0.96 

0.021 

-0.427 

0.395 

1.752 

0.92 


Where the equivalence ratio phi for H 2 /0 2 is given by: 


cHZ 

<j)o = 

2^cO° 2 


(30) 
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Chemical Kinetic Times For H 2 /0 2 

Step One Equations for H 2 /0 2 

The lean and rich step one H 2 /0 2 chemical kinetic time correlations are of the following 

form: 


r = A(P) a (<f> 0 ) b exp 


D 

T 


(31) 


where P is pressure in atm, <|) is the initial or overall equivalence ratio and T is the 
temperature in Kelvin. 

The initial or overall fuel oxygen ratio can be determined from either the initial or the 
mixture molar concentrations of cH 2 , cH 2 0 and c0 2 moles/cc. Assuming that the radical and 
intermediate concentrations are small, the initial cH 2 concentration is: 


cH 2° = cH 2 + cHIO 


(32) 


A summary of these correlations can be found in table 5. Parity plots for the lean step one 
H 2 /0 2 correlations can be found in figure 9. These parity plots show minimal scatter, which is 
consistent with the high R-squared values of the lean correlations. 


Table 5. — Step One H 2 /0 2 Chemical Kinetic Time Correlations 



A 

a 

D 

b 

Rsquared 

Lean 

2.632E-5 

-1.005 

24244. 

0.771 

0.975 

Rich 

1.397E-5 

-1 .024 

24799. 

1.141 

0.971 


(33) 

(34) 


Two Step H 2 /0 2 Chemical Kinetic Times 

Step Two Equations for H 2 /0 2 

The following form of equation was used for all lean and rich H 2 /0 2 step two 
correlations: 


t = A(P) a (cfuel) b (chlOy {(f) o) d exp 


E_ 

~T 


(35) 


Table 6 provides a summary of these correlations, and lean parity plots can be found in 
figure 10. 
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Table 6. — Step Two H2/O2 Chemical Kinetic Time Correlations 



A 

a 

E 

b 

c 

d 

Rsquared 

Lean 

1 .1 14E-13 

0.645 

8028 

- 1.000 

-0.664 

1.785 

0.358 

Rich 

2.196E-24 

2.544 

25388 

-3.079 

0.035 

4.24 

0.315 


(36) 

(37) 


Comparison of Kinetic Times for Hydrogen and 
NO x for Multiple Fuels 

The multiple fuel kinetic times at 10 atm and 2000 K are shown for step one in figure 1 1 and 
for step two in figure 12. The Jet-A and methane predictions were obtained from reference 2. For 
step one the kinetic times are one millisecond for lean mixtures except for Jet-A which was 
10 ps. The Jet-A mechanism consisted of a irreversible breakup to CH and C 2 H 2 followed by 
oxidation of the products. These times are the average startup times computed when the water 
concentration is less than 1 x 10 -20 . Without nitrogen, the rate for H 2 /O 2 is faster than FE/air. Step 
two, figure 12, is computed with an inefficiency of 0.5, or 50 percent reacted. After the radicals 
form, the rates are much faster being lps for H 2 /O 2 at an equivalence ratio of 0.6. The smaller the 
chemical kinetic time, the faster is the reaction. The chemical kinetic times for NO x step one and 
two (1NEFF = 0.5 see below) are shown in figure 13. The NO x production rate with hydrogen 
was two orders of magnitude slower than for Jet-A or methane. 


Computation of the Fuel and Water Concentrations for Step Two 

In order to use the step two instantaneous correlations, one needs to know the instantaneous 
fuel and water concentrations and the temperature and pressure. Remember, these correlations 
were based on the model concentrations for H 2 and H 2 O. The concentrations can be easily 
computed from the inefficiency, INEFF, where 

C H 2 

INEFF = (38) 

cH° 

2 

If cH° , the initial hydrogen concentration is unknown it equals: 


cH° = cH +cH O (39) 

222 v ' 

or can be computed from the initial f/a or initial equivalence ratio 4> 0 from: 

(40) where f/a=(|)o*0.0292 

(41) 


O _ f/a*29/2*p 
cHl ~ (1 + / /a *29/2) 

p = — moles / cc 

82.056 *T 
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Then, from equations (38) and (39), the concentrations are: 


cH 2 =cH° *(INEFF) 


cH 2 0= cH° -cH 2 = cH° *(1 -INEFF) 


cO = c cfi — ■■ cH iO 
2 ^2 2 


cN° = 

2 0.21 


0 79 * cO° 


(42) 

(43) 

(44) 

(45) 


The initial equivalence ratio (j) 0 for either H 2 /air or H 2 /0 2 is given by: 


(j)o — 


1 *cH° 2 

2 cO° 2 


(46) 


These values were substituted into step two kinetic time constants equations at 10 atm 
2000 K and an INEFF = 0.5 and plotted in figure 13. 

The step two chemical kinetic times can be computed from equations (24) and (35) above 
for a Tau given by the values in tables 3 and 6 at conditions of P = 10 atm, T = 2000 K, and 
INEFF = 0.5. 


Comparison of the Model With the Full Mechanism Kinetics 

The predictions of the full mechanism (H 2 , H,0,OH, and H 2 0) are plotted in figures 14 and 
15, and then our model predictions are superimposed on the calculations. The model predictions 
are slower than the full mechanism predictions. We have not included the intermediate species in 
our predictions, but we correlate the chemical kinetic times to the overall species concentrations. 
The inefficiency from the above equation (38) is integrated using the chemical kinetic time, Tau 

(<>• 


INEFF = 


cH 2 

o 

cH 2 



(47) 


where Taui = Taui at the start or 

Taui = Tau 2 when cH 2 0 > 1 x 10 _2 °. 

This single equation was integrated by backward differencing for a time step of 1 x 10 6 ms. 
Even at the first time step of 1 * 10” 6 sec, the water concentration is at 1 x 10 -13 , so only step two 
was used in these calculations. One could have used a Runge Kutta scheme, but this simple 
example could be quickly implemented on an Microsoft Excel spreadsheet. The chemical kinetic 
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time for this condition is about 10 jlis so the integration requires only about 3000 steps. The 
agreement is fair for the three species of interest. As expected the OH, H, and O increases at 
almost the same rate as the H 2 O, so it is a good indicator of the reaction rate. 


Computation of NO x PPM 

One can use the correlations for computing premixed NO x levels for well mixed combustors. 
The NO x in parts per million (ppm) or the emission index (El) can be computed using 
equation (6). 


NOx 


mole 

V cc J 


( 6 ) 


• NOx 


The chemical kinetic time tnox ms cc/mole was integrated for a stoichiometric equivalence 
ratio in figure 14. A NO x concentration of lxlO -9 moles/cc is predicted. 

The concentration in ppm is calculated as follows: 

^ moles \ 


NOx 


NOx( ppm) 


cc J 


moles 


x 10° where p 


v cc j 


RT 


6.09x10 


-5 


(48) 


The NO x concentration in ppm is then 16.42 for one millisecond residence time. 
This result is consistent with the data. 


Auto Ignition Times for H 2 /Air 

Hydrogen reacts explosively at low temperatures and pressures as shown by the dashed line 
in figure 16. Mixtures to the right of the curve explode (produce a rapid pressure rise for a 
constant volume reactor), while mixtures to the left are stable. This was computed by A.K. 
Oppenheim (ref. 8) using a full mechanism similar to the one used in this report with a one 
second reaction time. 

Auto ignition times using the simple model and a given formula were calculated and 
compared for the H 2 /air system. The auto ignition time for the simple model is based on the 
recommendation of reference 9, where the time required for ignition is for 5 percent of the fuel to 
react. 


Fuel = F o e 


(49) 
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(50) 


IfF “%, = °' 95 to = 0.051 3 and = 0.0513* r,„„. 


Note that r fuel decreases as f/a decreases. We used an equivalence ratio of 0.2 (fastest rate 
from figure 11) with the lean equation. 

We limited our correlations to pressures above 0.5 atm. The model above 0.5 atm fits the 
data well, but it will not capture the “S” shape of the classical data. At low pressures and 
temperatures, our model would overpredict the stability of the mixture. 

The auto ignition delay time was also compared to the correlations of reference 7, shown in 
figure 17. The formula for calculating the auto ignition time for H 2 (ref. 7) is given by: 


zPimilli sec onds ■ atm ) = 8 .E 6 exp 


9600 

T 


(51) 


Figure 17 shows a plot of the auto ignition time versus temperatures that the auto ignition 
time given by the simple model compares well with the data in reference 8 for x 2 and 
INEFF = 0.95 at the higher temperature. These correlations had a larger activation energy than 
the expressions of reference 8. 


Conclusions 

A simplified kinetic correlation scheme for H 2 /air and H 2 /O 2 fuels resulted in a two-time step 
correlation that calculates chemical kinetic times for fuel and NO x . These models are simple 
enough to use in spreadsheet calculations and they represent the data well. These chemical 
kinetic time equations can then be used in a numerical combustor code to compare the chemical 
kinetic time with the turbulent mixing time at every grid point and for every iteration step. 

Strong step one (R 2 > 0.9) and mild step two (R 2 ~ 0.4) H 2 correlations were developed (R 2 > 
0.9) for all relations, both the H 2 /air and H 2 /O 2 correlations. However, because we are trying to 
correlate so many values over a wide range of conditions, scatter error exists. These model 
equations are believed to be extremely useful in the comparison of kinetic reaction and turbulent 
mixing times and in the computation of kinetic rate results. In our model we have only two 
equations and four major species, therefore the step size can be large and the rates are no longer 
coupled. This should be very fast. 
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Appendix A 


Multiple Linear Regression With Microsoft Excel 


Performing Multiple Linear Regression on a Logarithmic Equation 


This regression technique may be used to develop a correlation between a dependent variable 
and one or more independent variables. First the equation to be used must be linearized. An 
example of an exponential equation used here is shown below. 


A = BC c D d exp 




J 


(Non-linear form) 


(Al) 


In (A) = In (B) + c ln(C) + d ln(Z)) + (Linear form) (A2) 

Columns of data containing the independent variables (natural log of C, natural log of D, 

1/T,) and the independent variable (natural log of A) were contained in an Excel spreadsheet. (It 
is easiest to have the independent variables adjacent to each other, followed by the dependent 
variable.) 

The multiple variable regression analysis is located in the Data Analysis Toolpak. The Data 
Analysis Toolpak must be added into the spreadsheet if it is not already running in Excel. In 
order to add it, select the ‘Add ins’ button from the Tools menu. Click on the Analysis Toolpak 
option and click OK to accept this choice. Then choose ‘Data Analysis’ from the Tools menu 
and double click on ‘regression’. Click on the ‘Input Y Range’ box and highlight the column that 
contains the logarithm of the dependent variable and press return. Click on the ‘Input X Range’ 
box and highlight the columns containing all of the independent variables. (In this case ln(C), 
ln(D) and 1/T). Press OK to begin the regression. The regression data will be contained in a new 
worksheet. The variable labeled ‘intercept’ will be equal to the natural log of coefficient B. The 
remaining coefficients (c,d, and e) will be given as X Variable 1,X Variable 2 and X Variable 3 
respectively. This process is quick and accurate for Excel 2002 and was used for all equations 
given in this report. Excel has to capability to handle one dependent variable and multiple 
independent variables. As many as 14,000 points were used in the regressions. Excel can regress 
up to 65000 points with one dependent variable and 8 independent variables. 
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Appendix B 


Modifications to GLSENS for Obtaining Output for the CH4 Computations 

C ALREADY DONE IN DIFFUN 
C MODIFICATIONS TO MAIN 

COMMON/marek/ foa, atauf , atauco, atauno, nc, eratio, edd, timi, nco 
C READ AND CONVERT INPUT, PERFORM PRE-KINETIC CALCULATIONS 
RXORDR= . TRUE . 
write (16,8) 

8 format ( ' fuel, nc, lnP,l/T,ln phi, ln(atau f ), atauf , millisec ' ) 
write ( 18 , 11 ) 

11 format ( ' nox, nc, lnP,l/T,ln phi ln(atau f ), atauf , millisec ' ) 
icount=0 
CALL KINP 
5 continue 
f o=0 . 

IF (NEXT) GO TO 55 
65 TIME1 = SECCPU (0.0) 

IF ( WSFLOW) WELSTR = .TRUE. 

C AFTER COMPLETION OF EACH TEST CASE, TABULATE THE AVERAGE CHEMICAL 
KINETIC TIMES 

pln=sngl (dlog (p) ) 

TLN=sngl (1 .dO/T) 
tl000=log ( 1000 . / sngl (T) ) 
f o=0 . 

phi=f oa/ 0.0292 
if (foa.gt.O.) FO=log (phi ) 

AF=0 . 

ACO=0 . 

ANO=0 . 

if (atauf . gt . 0 . ) AF=log (atauf) 
if (atauno . gt . 0 . ) ANO=log (ATAUNO) 
icount=icount+l 

write ( 0 , 31 ) i count , nc, nco, P, T, phi , timi 
C WRITE (14, 7) nc, PLN, TLN, FO, AF, ACO, ANO, timi 
WRITE ( 14 , 31 ) i count , nc, nco, P, T, phi, timi 
31 format (i7, 2il0, f 10 . 2, f 10 . 0, f 10 . 3, lp, lel3 .4) 

7 format (ilO, lp, 9el3 . 4 ) 
write (16, 7) nco, pin, tin, fo, af, atauf , timi 
write ( 18 , 7 ) nco, pin, tin, f o, ano, atauno, timi 
atauf=0 . 
nco=0 
nc=0 

atauno=0 . 

CALL RINP 
GO TO 5 

C MODIFICATIONS TO OUT2 OF LSENS COMPUTE STEP ONE AND STEP TWO FOR EACH 
PRINT STATION 

C h21ab . f 2/16/2005 
C COMPUTE MOLAR CONCENTRATIONS 
C ALREADY DONE IN DIFFUN 
420 f oa=eratio*0 . 02 92 


NAS A/TM— 2005-2 13619 


19 



Q O 


phi=eratio 
philn=log (phi) 
tn=timil 
trr=sngl (T) 
c write ( 6, 1605) time 

1605 format ( ' time ',el4.4) 

418 if (nc . gt . 0 ) goto 427 
write ( 10 , 1 ) 

1 format ( ' TITLE= P atm T K f/a coc h2 nox time') 
write (26, 444) dspnm ( 5 ) , dspnm ( 2 ) , dspnm ( 1 ) , dspnm ( 4 ) , dspnm ( 3 ) , dspnm ( 6 ) 
1 , dspnm ( 10 ) 

444 format ( ' timil ' , lp, 13 ( 3x, a8 , 2x) ) 

C calculate the initial conditions for the averaging 
nc=l 
nco=l 
NC2 = 1 
t0 = 0 . 
areaf=0 . 
tauf=0 . 
tauno=0 . 
areano=0 . 
atauf=0 . 
atauno=0 . 
do 600 i=l,ls 

if (dspnm(i) .eq. ' H2 ' ) ih2=i 
if (dspnm(i) .eq. 'H20 ' ) iih2o=i 
if (dspnm(i) .eq. 'NO ' ) ino=i 
if (dspnm(i) .eq. '02 ' ) io2=i 

600 continue 

STEP 1 INTEGRATION 

C set initial fuel concentration CH2MI 
CH2MI=SNGL (PRC (ih2) ) 

427 continue 

cnox=sngl (prc (ino) ) 
cfuel= (prc (ih2) ) 

CH20M=CH2MI-CFUEL 

if (ch2om. gt . 1 . d-20 . and. nco . gt . 8 ) goto 504 

IF (NCO/ 50*50 . EQ . NCO) write (10, 423) nco, P, T, phi, cfuel, h2osl 
ltimil 

423 format (14 , f 5 . 1 , f 8 . 1 , lp, 8el0.3) 

C calculate the initial conditions for the averaging 
if (nc . It . 2 ) goto 1500 
if (timil . gt . 1 .) goto 1500 
if (nc . eq . 1 ) t s tart = timil 
if (nc . eq . 1 ) t0=timil 

if (nc . eq . 1 ) tauf o=- (prc ( ih2 ) /W ( ih2 ) ) *1 . d3 
C if (nco . eq . 2 ) edd=15 . *tauf o 
if ( tn . eq . tO ) goto 424 

if (W ( ih2 ) . ne . 0 . ) tauf =- (prc ( ih2 ) /W ( ih2 ) ) *1 . d3 
if (nc.eq.q) atauf=tauf 
if (tauf . le . 0 .) goto 899 
nco=nco+l 

areaf=areaf+ (1 . /taufo+1 . /tauf) /2 . ^dble (tn-tO) 
timet=timxl.-t start 
atauf=dble (timet) /areaf 
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if (nco/l*l . eq . nco) write (ll,423)nc,P,T, phi , cf uel , 02 , tauf , atauf , 

ltimil 

tauf o=tauf 

xrate= (w ( ino) ) 

424 if (nco . eq. 1) tauni= (1 . DO/ (xrate) ) *1 . e3 
tauno=sngl ( 1 . DO/ (xrate) ) *1 . e3 
if ( tO . eq . tn) goto 339 
if ( tauno . it . 0 . ) go to 339 

areano=areano+ (1 . / tauno+1 . / tauni) /2 . *dble (tn-tO) 

atauno=dble (timet) /areano 

tauni=tauno 

c if ( (nc/10*10 ) . eq . nc) write (13,501)nc,P,T, cnox, tauno, 
c latauno , eratio, timil 
339 continue 

338 format (i7, lel2 . 3, 7el2 . 4) 
tO=tn 

504 continue 

431 format ( f 12 . 2 , 3el3 . 5 , f 8 . 3 , f 8 . 3 , el3 . 5 ) 

C do step 2 calculations 

IF (CH20M.LT. 1 .d-20) GOTO 899 

if (W(ih2) .ne. 0 . ) tauf=- (pro (ih2) /W (ih2) ) *1 . d3 

xrate= (w ( ino) ) 

tauni= ( 1 . DO/ (xrate) ) *1 . d3 

pl=dlog (p) 

tol=l . /T 

cfuel= (prc (ih2) ) 

h2o=sngl (prc (iih2o) ) 

h2osl=sngl (dlog (prc (iih2o) ) ) 

H20ML=dL0G (CH20M) 
fo=dlog (cfuel) 

if (W(ih2) .ne. 0 . ) tauf=- (prc (ih2) /W (ih2) )*X*d3 
if (tauf . le . 0 .) goto 899 
tal=dlog (tauf) 

if (tauni . gt . 0 . ) taln=dlog (tauni) 

NC2=NC2+1 

if (nc2/100*100 . ne . nc2 ) goto 899 
c if ( tauf . gt . 5 . . or . tauf . It . 1 . e-6) goto 899 
if (timil . gt . 1 . ) next= . true . 
if (timil . gt . 1 .) goto 1500 

write (10,423)nc2,P,T, phi , cfuel , h2o, ch2om 
1 , tauf, timil 

write (26, 447) timil, prc (5) , prc (2) , prc (1) , prc (4) , prc (3) , prc (6) , 
lprc ( 10 ) 

447 f ormat ( lp, 9el3 . 4 ) 

write (27 , 900 ) NC2 , timil , phi In, pi , tol , f o, h2oml , tal , tauf , ch2om 
write (28, 900 ) NC2 , timil , phi In, pi , tol , f o, h2oml , tain, tauni 
900 format (15, lelO . 3, lp, 9ell . 3) 
if (nc.gt.2)go to 502 
write ( 6 , 1602) ih2 , iih2o 
1602 formate ih2 ' ,i5,' iih2o ' ,i5) 
coci=prc (ih2) /ch2mi 
if (coci . gt . 0 . ) totau=-log (coci) 
write (6, 1600) timil, tauf , atauf , h2ol , h2oml 
1600 format ( ' timil T , e 1 2 . 3 , T tauf T , e 1 2 . 3 , T atauf f ,el2.3, 

1 1 h2ol ' , el2 . 3 , ' h2oml T ,el2.3) 
write ( 6 , 1603) ch2mi , cfuel , ch2om 
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1603 format ( ' ch2mi ',el2.3,' cfuel ' , el2 . 3 , 1 ch2om *,el2.3) 
write (6, 1601) coci, totau 

1601 format ( ' h2/h2i ',el2.3,' timil/tau overall ',el2.3) 
899 continue 

501 format (i4, f7 .4, f7 4el2 . 5 , f 6 . 3 , f 9 . 3 ) 

1500 continue 
nc=nc+l 

if (nc . gt . 3 ) goto 502 
DO 435 I J= 1 , MAX 

IF (IJ.GT.LS.OR.IJ.GT.LR) GO TO 435 
TCON (IJ) =SNGL (PRC (IJ) /W (IJ) ) 

FMOL=SNGL (SIGMA ( I J) *MIXMW) 

WRITE (LWRITE, 175 ) DSPNM(IJ) ,PRC(IJ) f FMOL f W(IJ) 

430 WRITE (LWRITE, 185) I J, RATE (IJ) , PRX (IJ) , EQUIL ( I J) 

C430 continue 
435 CONTINUE 

502 IF (WELSTR) GO TO 446 
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Appendix C 


Run Input for GLSENS 

Run input for H2 for GLSENS TYPICAL FOR STEP ONE AND TWO 
H2/AIR INPUT 
C02 AR 

TIME PRESSURE H2 

&prob TCON= . TRUE . , CT0=0 . 5 , 

IPRINT=1 , end=l . e-4 , send 

Sstart time=0 . ,ERATIO=0 .3, SCC=0 . 0, SCOX=0 . , SCH=2 .0, 

NOXRAT=3 .76, ARAT=0 . , CRAT=0 . , T=1 0 0 0 . , Send 
END 

& solver mxstep=20000 , emax=l . e-15, atolsp=l . e-15 , maxstp=l 90000 , Send 
I NIT 

ALLSP END 
SENSVAR 
ALLSP END 
REAC 
FINIS 


H202 INPUT FOR H2/02 STEP ONE AND TWO 
TIME PRESSURE H2 

&prob TCON= . TRUE . , CT0 = 0 . 5 , 

IPRINT=1 , END=1 . E-4 , &END 

&start time=0 . ,ERATIO=0 . 1, SCC=0 . 0, SCOX=0 . , SCH=2 .0, 

NOXRAT=0 . 0 , ARAT=0 . , CRAT=0 . , T=1000. , Send 
END 

& solver mxstep=20000 , emax=l .e-15, atolsp=l . e-15 , maxstp=l 90000 , &end 
INIT 

ALLSP END 
SENSVAR 
ALLSP END 
REAC 
FINIS 
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Figure 1.— Magnussen mixing model. 



Figure 2. — H 2 /air T4 parity (lean) R 2 = 0.936. 
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Figure 3.— N0 X mole fraction equilibrium parity (lean) R 2 = 0.977. 
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Figure 7. — Kinetic H 2 /air N0 X parity plot step two (lean) R 2 = 0.976, 
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Figure 9,— Kinetic H 2 /O 2 parity plot step one (lean) R 2 = 0.956. 
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Figure 10.— Kinetic H 2 /O 2 parity plot step two (lean) R 2 = 0494. 
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Figure 1 1. — Kinetic i fuel at 10 atm and 2000 K for fuel reactions step one (lean and rich), time = 0. 
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Figure 12, — Kinetic t H 2 at 10 atm and 2000 K step two (lean and rich), INEFF = 0,5. 
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Figure 1 3. — Kinetic x N0 X at 10 atm and 2000 K for NO x reaction step one multiple fuels, time - G f x 2 for 
H 2 /air (lean and rich), and INEFF = 0.5. 
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Figure 16. — Explosion limits for H 2 at low pressure. 
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